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DEFINITION OF SYMBOLS 


SYMBOL 

A 



c 


D 


DEFINITION 

Aircraft reference wing area (m^) 

Drag coefficient, D/(l/2)pV^A (dimensionless) 

Drag coefficient at zero angle of attack (dimensionless) 

Drag derivative (rad ^) 

_2 

Drag derivative (rad ) 

Drag derivative (rad”^) 

Lift coefficient L/(l/2) pV^A (dimensionless) 

Lift coefficient at zero angle of attack (dimensionless) 
Lift derivative (rad ^) 

Lift derivative (rad ^) 

Lift derivative (rad ^) 

Lift derivative (rad ) 

Pitching-moment coefficient M/(l/2)pV^Ac (dimensionless) 

Pitching-moment coefficient at zero angle of attack 
(dimensionless) 

Pitching-moment derivative (rad“^) 

Pitching-moment derivative (rad~^) 

Pitching-moment derivative (rad ^) 

Pitching-moment derivative (rad”^) 

Wing mean aerodynamic chord (m) 

Aerodynamic drag force (N) 


vi i i 



SYMBOL 


DEFINITION 


F Net propulsive thrust (N) 

g Acceleration due to gravity (m/s^) 

GS+ Flight path deterioration parameter for ascent above glide 

slope (dimensionless) 

GS- Flight path deterioration parameter for descent below glide 

slope (dimensionless) 

lyy Aircraft moment of inertia about lateral y axis (kg - m^) 

L Aerodynamic lift force (N) 

m Aircraft mass (kg) 

M Pitching moment (kg m^/s^) 

q Pitch angular velocity (rad/sec) 

t^ Engine characteristic spool -up time (sec) 

V Magnitude of airspeed (m/s) 

Magnitude of initial airspeed (m/s) 

V+ Flight path deterioration parameter for airspeed deviation 

above (m/s) 

V- Flight path deterioration parameter for airspeed deviation 

below (m/s) 

W Magnitude of horizontal wind component (head wind +, tail 

wind -) (m/s) 

W- Magnitude of vertical wind component (down draft +, updraft -) 

(m/s) 

X Horizontal coordinate (m) 

X Aircraft ground speed (m/s) 

z Vertical coordinate (m) 

z Aircraft vertical velocity (m/s) 
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SYMBOL 


DEFINITION 



Elevator deflection angle (rad) 

A 

Denotes differential 

0 

Aircraft pitch angle (rad) 

P 

Air density (kg/m^) 

"ph 

Other 

Phugoid frequency (sec"^) 

(■) 

First time derivative 

(••) 

Second time derivative 

( )f( 

Nominal (trim) value 

( )o 

Initial value 


X 



CHAPTER 1. INTRODUCTION 


Several aircraft accidents and near accidents have been officially 
attributed to encounter with low-level wind shear during the approach 
and takeoff phases of operation. The most notable accident with wind 
shear as a major cause involved the Eastern 66 Boeing 727 (B727) which 
crashed on approach in a thunderstorm wind environment at John F. 

Kennedy (JFK) International Airport on June 24, 1975, resulting in 113 
fatalities [1].* 

Wind shear is a spatial or temporal gradient in wind speed and/or 
direction. Strong shears are generally associated with the presence of 
cold and warm fronts and thunderstorm cells. An analysis of the Eastern 
66 accident by Fujita and Caracena [2] revealed that Eastern 66 was 
forced above and then subsequently below the glide slope when it flew 
directly into a downburst cell consisting of a strong down draft and a 
head wind shearing to a tail wind. 

While Fujita blamed the accident on the presence of strong down 
drafts, several subsequent studies by McCarthy, Blick, and Bensch [3], 
Frost and Crosby [4], Frost and Reddy [5], and Turkel and Frost [6] 
focused mainly on fixed-controlled and pilot-controlled aircraft response 
to longitudinal wind shears, which excite the aircraft's phugoidal (long 
period) oscillatory mode. The phugoidal time period, which depends on 
the aircraft's airspeed (approximately 70 m/s for a B727 on approach), 
is about 38 sec. At an approach ground speed of 70 m/s, an aircraft 
flies 2.7 km in 38 sec. This is well within a thunderstorm downburst 
cell's horizontal extent [2]. It is therefore possible for an aircraft 
on approach to encounter a thunderstorm downburst cell and experience a 
wave-like head wind to tail wind shear at approximately the aircraft's 
phugoidal frequency. During the thunderstorm at JFK, several aircraft 


*Numbers in brackets correspond to similarly numbered references in 
the list of references. 
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landed with varying degrees of success. Seven aircraft experienced 
little or no problem; one had extreme difficulty; one executed a missed 
approach; and Eastern 66 crashed. Fujita theorized that these aircraft 
passed through various sections of several downburst centers crossing 
the approach path at various times and that Eastern 66 flew directly 
into the center of a strong cell. 

Since the accident, several wind shear detection and warning systems 
and operational procedures have been proposed to prevent further acci- 
dents and alert flight crews and air traffic controllers to the presence 
of wind shear. Examples include LLWSAS (low level wind shear alert 
system) [7], PODS (pressure jump detector system) [8], and the airspeed- 
ground speed operational procedure [7]. However, none of the proposed 
systems provides a warning of severe wind shear located along the 
approach path prior to encounter by the aircraft. The LLWSAS and PJDS 
provide low level warning at the ground at the airport and in its sur- 
rounding areas. The airspeed-ground speed monitoring procedure requires 
the aircraft to be already in the wind shear and to have ground speed 
information readily available. 

A wind shear warning and detection system capable of providing ade- 
quate warning of severe wind shear has been proposed* by McCarthy [9]. 

The system would consist of a Doppler radar, a microprocessor, a data 
link, and displays. The Doppler radar (Figure 1) would be located 
near the runway, aimed along the approach or takeoff path, and would 
be programmed to rapidly measure horizontal wind components at 150 m 
slant length intervals between the ground and an altitude of 500 m. The 
measured flight path longitudinal wind field and the aircraft configura- 
tion designation would be entered into a microprocessor which would then 
integrate the equations of motion to predict the aircraft's probable 
deviation from its desired flight path and to compute a set of parameters 
which describe the aircraft's deviation from its desired flight path. 
These so-called flight path deterioration parameters, which will be 
discussed in detail below, are used to quantify the aircraft's trajectory 
behavior as a single set of numbers. The flight path deterioration 
parameters would then be transmitted via the data link to the air traffic 
controller and pilot's data display (Figure 2). 
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Figure 2 Doppler radar wind shear warning and detection system. 




The Doppler radar is an excellent device for wind speed detection 
since it can measure the longitudinal wind for long distances along 
flight paths. One Doppler beamed along the glide slope only measures 
the wind velocity component that is aligned with the Doppler beam, which 
is approximately equal to the horizontal wind component along the flight 
path. Because the vertical velocities along the glide slope are diffi- 
cult to measure with Doppler radar and because the longitudinal winds 
directly affect the aircraft's longitudinal stability, attention is 
focused on using only longitudinal wind data as input to warning system 
software. The vertical velocity components should not be ignored, 
however, since they clearly influence the longitudinal wind field. 

The aim of this study is: 

1. To utilize an aircraft trajectory program to examine aircraft/ 
pilot response through wind shears including longitudinal sine 
waves, S-shape waves, 1 - cosine vertical winds, and combina- 
tions at various frequencies and amplitudes as approximations 
to the winds encountered in a thunderstorm downburst cell. 

2. To determine if the control system algorithm and aircraft 
trajectory program combination gives an accurate representa- 
tion of the behavior of the real pilot in a flight simulator 
and hence real aircraft when subjected to the input wind field 
models described above. 

This report will discuss the trajectory program, the pilot/control 
system models, the quantitative flight path deterioration parameters, 
the trajectory program results, the flight simulator results, and will 
present conclusions and recommendations. 



CHAPTER 2. AIRCRAFT TRAJECTORY PROGRAM 


2.1 Program Specifics 

The aircraft trajectory program is a two-dimensional, three-degree- 
of-freedom model with translational motion in the longitudinal (x) and 
vertical ( 2 ) planes and rotational motion about the lateral (y) axis. 

A two-dimensional model is sufficient for a warning and detection system 
since we are looking primarily at the effects of longitudinal winds along 
the glide slope which cause the aircraft deviations in the x-z plane of 
the glide slope. The effects of lateral winds on the aircraft's hori- 
zontal deviation of the glide slope have not been included in the air- 
craft trajectory model since the lateral wind component cannot be 
measured by the single Doppler radar. The three-degree-of-freedom 
warning system program developed in this study requires little computer 
storage capacity and has execution times compatible with real-time field 
operations. It should be noted that the aircraft trajectory program was 
derived from the previous model of Turkel and Frost [6]. The program 
was compacted and the pilot/control system was both improved and simpli- 
fied. The new model has a run time less than 3 sec on a VAX 11/780 as 
compared to 50 sec for the previous model. The equations of motion are 
discussed in Section 2.2. 

The aircraft trajectory program is adaptable to any aircraft. The 
aerodynamic lift, drag, and moment coefficients and other data (i.e., 
wing area, weight) and the initial control settings must be input by the 
user. In a real-time system, these data may be stored by aircraft type 
to be looked up for a specific aircraft entering the warning zone. In 
addition, the program can be run in the takeoff or landing configuration. 
The Doppler-measured longitudinal winds can be entered as a function of 
horizontal distance from the* runway threshold in a subroutine to be 
utilized in the simulation of the aircraft's trajectory. After the simu- 
lated trajectory through the Doppler-monitored warning zone is completed, 
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the flight path deterioration parameters based on airspeed and flight 
path deviation can be calculated and any potential hazards due to the 
longitudinal wind field can be identified by the air traffic controller. 
These parameters are described in Chapter 3. 

The trajectory program control systems can simulate aircraft control 
by autopilots or pilots of varying degrees of skill. Engine spool -up 
lag and pilot thrust command delay are included. These thrust command 
lags due to the pilot response delay and the turbomachinery response 
characteristics are important factors in a longitudinal wind shear 
encounter where the aircraft needs a large power increase in a short time 
to stay out of danger. The control systems are described in detail in 
Section 2.3. 


2.2 Equations of Motion 

The forces acting on the aircraft and angular relationships for 
the two-dimensional , three-degree-of -freedom model are presented in 
Figure 3. The equations of motion in earth coordinates are: 

1) Acceleration along longitudinal axis 

"F L./\D /\ 

X = — cos e - — sin(Y) - — cos(y) 
m m m ' ' 

2) Acceleration along vertical axis 

z = ^ sin e + ^ cos(y) - ^ sin(Y) - g 
m m ' ^ m \ i ^ 

3) Acceleration about lateral axis 
M 


e = 


I 


yy 


where aircraft lift, drag, and moment are given respectively by 
L = 1/2 pV^C^A 


D = 1/2 pV^CpA 
M = 1/2 pV^C^Ac 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 
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Figure 3 Relationship between the various forces acting on an aircraft. 


The magnitude of the relative wind vector, i.e., the airspeed, is given 
by 

V = [(i + W^)2 + a + (7) 

where x and z are the aircraft ground speed and vertical velocity rela- 
tive to the earth. W is the longitudinal wind speed component and is 

A 

defined as positive when opposing the aircraft and is hence called the 
head wind. Similarly, is the vertical wind speed component and is 
defined as positive in the negative z direction and is commonly called 
the down draft. 


= tan 


The angle y is given 
-1 


+ W^l 


(8) 


and angle of attack is given by 

a = 0 + (- 7 ) = 0 - Y (9) 

The aerodynamic lift, drag, and moment coefficients in Equations 4 
through 6 are given by 


Cj^ - C|^ + C|^ (ct) + C|^ 2V 

‘"0 '~a q a 


( 10 ) 
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(11) 


Cp = Cp + Cp (a) + Cp y) * Cp (6^) 

0 a 

= ‘^n, ^ ^ fv ^ Iv (^2) 

0 ct o^ q a 

The aerodynamic data used in this study represents a B727-type air- 

craft in the landing or takeoff configuration with the flaps deflected 
30° and the landing gear extended. This data is given in Appendix A. 

The unsteady terms for q and a have been included because the 
aerodynamic forces on the aircraft take some finite amount of time to 
adjust to the changing wind conditions and aircraft configurations which 
occur when an aircraft encounters a thunderstorm wind field. The effect 
of these terms on aircraft dynamics is discussed in Etkin [10]. 

2.3 Pilot/Control System Models 

The logic for three different control system models has been devel- 
oped for simplified and meaningful simulation of controlled aircraft 
flight. The three control systems model the behavior of a typical pilot 
and two different autopilots. All three systems utilize throttle com- 
mand rate to control thrust and thrust rate. The elevators of the B727 
are used as control surfaces for pitch, vertical velocity, and flight 
path angle control. All three models use fixed gains for thrust rate 
and elevator angle adjustments during flight in varying winds. These 
gains can be adjusted for simulation of pilots of varying skill or 
autopilot control. Engine spool-up lag time due to turbomachinery lag 
and pilot thrust command delay have been included in the control system 
since instantaneous thrust command by the pilot does not provide instan- 
taneous desired thrust. The pilot initiates a throttle rate command to 
control thrust and thrust rate. Because of the turbomachinery and pilot 
response lags, a predictive control scheme has been included. The 
purpose of this scheme is to provide some degree of lead since even the 
most responsive pilot must provide some lead in his throttle conmand 
because of turbomachinery lag. The turbomachinery lag time for an 
engine -to reach full thrust from normal approach rpm can be as small as 
2 sec and as high as 6 sec from idle rpm [11]. No characteristic lag 
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time was applied to elevator angle control since elevator response times 
are usually very small. 


The B727 has three engines which deliver a total maximum thrust of 
187,000 N (42,000 lbs) at takeoff at sea level and a minimum of 13,350 N 
(3,000 lbs) at forward idle. The maximum rate of thrust increase or 
decrease is on the order of 36,000 N/s (8,000 Ibs/s) [12]. The thrust 
time constant (t^) which describes the nonlinear time for thrust to 
advance from low power to full power at full setting could be as small as 
2 sec [11]. 

The thrust control lag has been modeled assuming that the difference 
between the change of thrust and change of throttle setting can be 
expressed as a decaying exponential of the form 


dF = 



(13) 


At time t = t, before the throttle is moved, the change in thrust is 
zero. At the time when the throttle is moved, a thrust change command 
(dAF) is produced. At time t = t + t^, dF produced by the throttle move- 
ment input is (1 - 1/e) of the desired dAF. At t = t + 2t^, the change of 
thrust dF is (1 - l/e^) of the desired dAF and exponentially approaches 
dF as t ^ “. 


The equation for thrust from t = 0 to t = t allowing for continuous 
throttle movements included between t = 0 and t = t is given by 




- exp 



(14) 


where dAF/dx is the rate of change of throttle setting. 

A more convenient form for Equation 14 than the integral form is a 
differential equation. Therefore, differentiating Equation 14 with 
respect to x yields the following expression for F (rate change of 
thrust) 



( 15 ) 


10 



so that 


t 

F + = F„ + f ^ <Jt 

C 0 J ClT 

0 

Differentiating Equation 16 further yields 


F + V = ^ (t) 


which can be written 


( 16 ) 


(17) 


F = 


dAF 

dr 


(t) - f 


(18) 


J A 17 

where (t) is the rate change of throttle- setting resulting from 
throttle movement by the pilot. 


Three different control system models were investigated to simulate 
pilot or autopilot control. The first control system, model I, uses 
throttle setting control to initiate thrust rate control for controlling 
changes in aircraft ground speed. Elevator control is applied to correct 
for changes in vertical velocity and pitch angle. The control equation 
for thrust rate is given by 


u = ^ [(x + t x) - Xn|] 

3X 


(19) 


d A F 

where u = (t) is the rate of change of throttle setting commanded by 

the pilot's change in throttle position. In Equation 19, Uj^ corresponds 
to the trim throttle position change command which is always equal to 
zero at trim. The partial derivative, 9u/9x, is a fixed gain for pilot 
throttle rate command (by throttle position change) for changes in 
predicted aircraft ground speed (x + t^x) from the trim condition ground 
speed X|^. The utilization of predicted ground speed enables the pilot 
or autopilot to compensate for the turbomachinery delay and for pilot/ 
autopilot response lag. In this case, t^ = 5.0 sec (a relatively slow 
turbomachinery/pi lot response lag time constant). The 9u/9x gain used 
in this control system logic was set at -1,000 (N/s)/(m/s) which commands 
a ^j^jl,000 N/s thrust rate by throttle position movement for a 1 m/s 
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increase (decrease) in ground speed. It should be noted that u = 0 when 
the throttle is at idle or full -forward settings. 

The control equation for elevator setting is given by 

36p 96p 

\ ^ ^ ( 20 ) 

where 6r. is the desired elevator setting and 6c is the trim elevator 
setting. The fixed gains 86^/az and 96^/96 correspond to elevator angle 
commands due to changes in vertical velocity and pitch angle from the trim 
values, respectively. The 96^/9z gain used in this control system model 
was set at 0.01 rad/ (m/s) corresponding to an elevator angle command to 
pitch nose up (nose down) with 0.01 rad of elevator angle with 1 m/s 
increase (decrease) in sink rate. The 96^/99 gain used in this control 
system model was set at 0.1 rad/rad corresponding to an elevator angle 
command to pitch nose up (nose down) with 0.1 rad of elevator angle with 
1 rad of nose pitch down (nose pitch up). It should be noted that the 
pilot/autopilot response characteristics can be altered or "fine-tuned" 
by changing the magnitudes of t^ or any of the three gains. Setting the 
three gains to zero corresponds to a fixed control mode of operation. 

The second control system, model II, uses throttle setting control to 
initiate thrust rate control for controlling changes in the aircraft 
airspeed. Elevator control is applied to correct for changes in flight 
path angle and pitch angle. The control equation for thrust rate is 
given by 

U = [(V + t^^) - V,^] (21) 

where u = (t) is the rate change of throttle setting commanded by 

the pilot by change in throttle position, and Uj^ = 0 as noted previously 
in Equation 19. The partial derivative, 9u/9V, is a fixed gain for 
pilot throttle rate command (by throttle position change) for changes in 
predicted aircraft airspeed (V + t^^) from the trim condition airspeed 
V|^. The 9U/9V gain used in this control system logic was set at -2,000 
(N/s)/(m/s) which commands a ^^^^2,000 N/s thrust rate by throttle 
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position movement for a 1 m/s increase (decrease) in airspeed and t^ = 
2.0 sec. 

The control equation for elevator setting is given by 


96p 96p 


( 22 ) 


where 6r.» 6r » 96r-/9z, 96^/90 are defined as previously in Equation 20. 

t t|^ t t 

Flight path angle is maintained by vertical velocity control through 
elevator angle control. To maintain a desired flight path angle of 
approximately 3 degrees, the pilot must attain a desired vertical 
velocity equal to: 

= -x(0.052) (23) 


where x is the instantaneous aircraft ground speed and 0.052 is approxi- 
mately tan (3°). The gains 96p/9z and 96^/96 were set to 0.01 rad/(m/s) 
and 0.1 rad/rad, respectively. This system works well when the aircraft 
is held relatively close to the glide slope beam and the flight path 
angles are kept within reasonable limits. However, when the aircraft 
deviates a large distance above or below the glide slope for a relatively 
long time period, the logic of tracking a 3° glide slope beam is insuf- 
ficient for landing at the touchdown point. 

The third control system, model III, controls airspeed by throttle 
rate application, with touchdown point and pitch attitude controlled by 
elevator angle control. The control equation for the thrust rate is 
given by Equation 21 of the previously mentioned airspeed control system. 
A turbomachinery/pilot response lag time constant t^ was set at 3 seconds 
and the thrust rate gain 9u/9V was set at -1,000 (N/s)/(m/s). The 
control equation for elevator angle control is the same as that in the 
previous control system given by Equation 22. However, in this case, 

Z|^, the desired vertical velocity is determined by the aircraft's posi- 
tion above or below the glide slope. If the aircraft is above the glide 
slope, Zp is determined by the flight path angle needed to land at the 
horizontal touchdown point Xj^. This relationship is given by: 
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where x is the aircraft ground speed and z and x are aircraft vertical 
position and horizontal position, respectively. If the aircraft falls 
below the glide slope, the control system schedules a 1° ascent flight 
path angle to reintercept the glide slope, rather than fly a shallow 
flight path angle to touchdown at Xj^. The desired ascent rate for a 1° 
glide slope intercept is given by 

Zp = x(0.0174) (25) 

where x is aircraft ground speed and 0.0174 is approximately tan (1°). 

The elevator control gains 36^/9z and 96^/96 were set at 0.015 (rad)/(m/s) 
and 0.1 rad/rad, respectively. 
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CHAPTER 3. FLIGHT PATH DETERIORATION PARAMETERS 


In order to assess the potential severity of a longitudinal wind 
shear hazard existing along the flight path, some quantitative parameters 
are needed to describe the response of the aircraft/pilot system. These 
parameters are called the flight path deterioration parameters (FPDP) and 
they describe the aircraft trajectory by glide slope deviations and air- 
speed deviations from reference airspeed. The Doppler-measured longitu- 
dinal glide slope winds are fed into the microprocessor containing the 
aircraft trajectory program and a simulated trajectory is calculated. 

The FPDP would then be computed by the program to describe the aircraft's 
deviations from the glide slope and reference airspeed. These parameters 
would be displayed to the controller or pilot and a decision as to 
whether or not to attempt takeoff or landing along the glide slope moni- 
tored by the wind shear warning system could be made before the aircraft 
encounters the monitored area. Naturally, there exists a need to deter- 
mine the safety limits of these parameters so that the pilot or con- 
troller can make a go/no-go decision based on the parameters displayed by 
the warning and detection system. It should be noted that the parameters 
used in this study are basically for analysis purposes--they may not be 
the most informative data that a pilot or controller can relate to a 
safe approach or takeoff. However, it is anticipated that the final 
parameters used in the actual system will be based on similar analytical 
parameters. 

The proposed flight path deterioration parameters for flight path 
and airspeed deviation were devised using the following logic. Normal- 
ized altitude, HP/HG, was chosen as the flight path deviation parameter 
where HP is the height of the aircraft above the ground and HG is the 
height of the glide slope above the ground. This formulation was chosen 
over aH = HP - HG because it accounts for the fact that altitude devia- 
tions become more serious at lower altitudes than at higher altitudes. 
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The airspeed deviation parameter chosen was simply the deviation of air- 
speed from the reference value. However, in both flight path and 
airspeed parameters, the positive deviations are examined separately 
from the negative deviations to avoid the cancellation errors. Although 
an averaging technique is used, it is believed that this method provides 
more meaningful results than root mean square averaging. 

In summary, the FPDP used in this analysis and illustrated in 
Figure 4 are: 


1. Flight path deviation 

a) GS+ = ^ I 

" 1=1 '^‘^1 


HP^ = height of aircraft above ground at point i 
HG^ = height of glide slope above ground at point i 

where HP^./HG^. above or on glide slope >1.0 

n = number of points aircraft is above glide slope 

The points are taken at 1 sec intervals. 

. m HP. 

b) GS- = i .1 hGT 
1=1 1 


where HP^/HG^ below glide slope < 1.0 
m = number of points aircraft is below glide slope. 



Figure 4 Concept of flight path deterioration parameters (FPDP). 
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2. Airspeed deviation 

a) v+ = i (V,. - 

V^. = aircraft airspeed at point i 
= reference airspeed 
where (V- - > 0 

j = number of points aircraft airspeed (V^. ) is above 

b) V- = 1 (V,. - 

where (V- - < 0 

k = number of points aircraft airspeed (V^- ) is below 

A detailed example illustrating the calculation of the FPDP is 
given in Appendix B. 
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CHAPTER 4. RESULTS OF COMPUTER AND SIMULATOR STUDIES 
4.1 Description of Test Plan 

In an effort to determine if the aircraft trajectory model simulates 
a real aircraft/pilot system, the results of the trajectory program were 
compared with a series of runs that were carried out in the B727 simula- 
tor at the NASA Ames Research Center (NASA Ames). The simulated aircraft 
was flown through three different wave-form wind models characteristic of 
a thunderstorm downburst cell environment (Figure 5). The simulator runs 
were designed to test aircraft/pilot response to longitudinal and 
vertical wind waves of varying amplitudes and frequencies. The test plan 
appears in Table 1 and a summary of the computer and simulator results 
that will be discussed below is included in Tables 2 and 3, respectively. 

4.2 Illustration of Aircraft Sensitivity to Longitudinal 
Wind Shears at Various Frequencies 

Before discussion of the results of the aforementioned computer and 
simulation tests, it is desirable to examine the response of the fixed- 
control aircraft (no pilot input) to winds described as longitudinal sine 
waves of constant amplitude but different frequencies using the aircraft 
trajectory program as developed in this study. This will serve to illus- 
trate the hypotheses of McCarthy, Blick, and Bensch [3], and Frost and 
Crosby [4] stating that fixed-control aircraft are most sensitive to 
longitudinal wind shears occurring at approximately the aircraft's 
phugoid frequency. The following tests were made using the aircraft 
trajectory model with an airplane characteristic of a B727. The aircraft 
was trimmed (flaps 30°) for flight along a 2.9° flight path with an 
approach airspeed of 75.1 m/s and an angle of attack of 4.7°. In all 
cases the wave is encountered at x = 0.0 m. Figure 6a shows the trajec- 
tory of an aircraft characteristic of a B727 encountering a 5.15 m/s 
(10 kt) full head wind to tail wind sine wave at the phugoid frequency 
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Head Wind 


Tail Wind 



Down Draft 


Updraft 



Figure 5 Wind models used to simulate a thunderstorm downburst cell. 
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TABLE 1 Test Plan for Simulator and Computer Runs. 



Ampl itude 

Frequency 

Aircraft trimmed for: 3.0° glide slope 



70.0 m/s airspeed 



63,958 kg {140,000 lbs) 



gear down, flaps 30° 



S-Shape head wind to tall wind shear wave 

5.15 tn/s (10 kts) 



10.30 m/s (20 kts) 

> 


15.45 m/s (30 kts) 

Up^ (=38 sec) 


15.45 m/s 

2up^ (=19 sec) 


15.45 m/s 

1/2 uip^ (=76 sec) 

Full sine wave head wind to tall wind shear 

5,15 m/s 

"ph 


10.30 m/s 

"ph 


15.45 m/s 

> 


20.60 m/s 

■"ph 


20.60 m/s 

^ “ph 


20.60 m/s 

“ph 

1 - cosine down draft 

5.15 m/s (17 ft/s) 

“ph 


10.30 m/s (34 ft/s) 

■^ph 


15.45 m/s (51 ft/s) 

“ph 


15.45 m/s (51 ft/s) 

2 “ph 


15.45 m/s (51 ft/s) 

V2 > 

Combinations: 



fS-shape, 0-5.15 m/s tall wind shear at 



ll - cosine, 5.15 m/s (17 ft/s) downburst at Upjj 



rs-shape, 0-10.3 m/s tall wind shear at Upj, 



ll - cosine, 10.3 m/s (34 ft/s) downburst at Up^ 



rS-shape, 0-20.6 m/s tall wind shear at Wp^ 



- cosine, 20.6 m/s (68 ft/s) downburst at uipp 



fS-shape, 0-20.6 m/s tall wind shear at 2 cupj^ 



ll - cosine, 20.6 m/s (68 ft/sl downburst at 2 «p^ 
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TABLE 2 Flight Trajectory Computer Run Summary 


Figure 


Wave 

Shape 

Direction 

Amplitude 

(m/s) 

Frequency 

^air 
(m/s) 
(min, max) 

a 

(Deg) 
(min, max) 

V+ 

(m/s) 

V- 

(m/s) 

i GS+ 

1 GS- 

1 '^/'^max* 

(min, max) 

6a 


Sine 

Longitudinal 

5.15 

1.0 

65, 88 

3.7, 6.1 


-5.93 

1.32 

0.56 

0.35, 0.35 

6b 

Fixed 

Sine 

Longitudinal 

5.15 

2.0 

67, 80 

4.2, 5.9 

2.78 

-3.68 

1.11 

0.76 

0.35, 0.35 

Be 

Fixed 

Sine 

Longitudinal 

5.15 

0.5 

71, 79 

4.4, 5.2 

1.92 

-2.37 

2.48 

0.13 

0.35, 0.35 

7a 

Pilot 

Sine 

Longitudinal 

5.15 


73, 79 

3.8, 5.4 

1.54 

-1.52 


0.84 

0.21, 0.47 

7b 

Pilot 

Sine 

Longitudinal 

5.15 


70, 79 

3.8, 5.9 

0,97 

-1.24 


0.97 


7c 

Pilot 

Sine 

Longitudinal 

5.15 

0.5 

71, 78 

3.8, 5.8 


-2.56 

1.30 


0.16, 0.51 

8a 

Pilot 

Sine 

Longitudinal 

5.15 

1.0 

63.7, 81.2 

3.4, 8.5 

4.01 

-3.57 

1.09 

0.78 

1, 0.83 

8b 

Pilot 

Sine 

Longitudinal 

10.30 


59.8, 89.6 

1.8, 10.0 


-4.78 

1.26 

0.83 

I, 1.0 

8c 

Pilot 

Si ne 

Longitudinal 

15.45 


56.9, 93.8 

1.1, 11.4' 

8.99 

-6.09 

1.40 

0.79 

I, 1.0 

10a 

Pilot 

Sine 

Longitudinal 

20.60 


54.4, 97.0 

0.6, 12.6 



1.50 

0.73 

I. 1.0 

10b 

Pilot 

Sine 

Longitudinal 

20.60 


45.4, 88.7 

0.1, Stall 

ill 



0.69 

I, 1.0 

10c 

Pilot 

Sine 

Longitudinal 

20.60 

0.5 

63.0, 84.0 

2.9, 8.7 

4.97 

-3.71 

1.77 

0.94 

I, 0.95 

12a 

Pilot 

1- Cosine 

Downburst 

5.15 

1.0 

67.0, 74.0 

5.0. 7.4 

1.97 

-1.63 


0.81 


12b 

Pilot 

1- Cosine 

Downburst 

10.30 

1.0 

69.3, 78.7 

3.8, 8.0 


-2.77 


0.79 


12c 

Pilot 

1-Cosine 

Downburst 

15.45 

1.0 

63.0, 80.1 

3.4, 8.7 

5.28 

-3.66 



I, 1.0 

14b 

Pilot 

1- Cosine 

Downburst 

15.45 

2.0 

60.5, 81.8 

3.0, 10.3 

5.24 

-3.89 


0.75 

I. 1.0 

14c 

Pilot 

1-Cosine 

Downburst 

15.45 

0.5 

66.6, 70.1 

6.2, 7.3 


-1,89 


0.69 

0.33, 1.0 

16a 

Pilot 

S-Shape 

Longitudinal 

5.15 

1.0 

68.3, 72.1 

5.6, 6.9 

1.13 



0.98 


16b 

Pilot 

S-Shape 

Longitudinal 



66.3, 74.4 

5.0, 7.6 

2.21 

-1.84 


0.95 


16c 

Pilot 

S-Shape 

Longitudinal 

15.45 


65.3, 76.4 

4.4, 8.0 

3.36 

-2.15 

1.24 

0.95 

I, 0.87 

18b 

Pilot 

S-Shape 

Longitudinal 

15.45 


62.1, 81.8 

3.1, 9.0 

5.63 

-3.79 

1.12 

0.83 


18c 

Pilot 

S-Shape 

Longitudinal 

15.45 

0.5 

69.2, 71.4 

5.8, 6.5 

0.85 


1.60 

0.98 
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Pilot 

1- Cosine 

Downburst 

5.15 

1.0 










5 -Shape 

Longitudinal 

5.15 

1.0 

66.2, 75.4 

4.6, 7.6 

2.53 

-1.93 




22a 

Pilot 

1 -Cosine 

Downburst 

10.30 

1.0 










S-Shape 

Longitudinal 

10.30 

1.0 

64.5, 80.5 

3.5, 8.2 

5.28 




I, 1.0 

22b 

Pilot 

1-Cosine 

Downburst 

20.60 

1.0 










S-Shape 

Longitudinal 

20.60 


54.3, 70.1 

6.0, 13.3 




0.67 


22c 

Pilot 

1-Cosine 

Downburst 

20.60 

2.0 










S-Shape 

Longitudinal 

20.60 

2.0 

49.8, 75.9 

6.3, 16.0 

2.89 

-10.39 


0.62 
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Auto- 

Pilot 

III 

Sine 

Longitudinal 


1.0 

56.9, 88.5 

-0.3, 12.7 

7.40 

-5.52 

1.44 

0.86 

I, 1.0 


*I = idle thrust. 






TABLE 3 


Figure 

Control 

System 

Wave Shape 

Direction 


9a 

Pilot 

Sine 

Longitudinal 

5.15 

9b 



II 

10.30 

9c 



II 

15.45 

11a 



II 


11b 


" 


20.60 

11c 





13a 

Pilot 

1-Cosine 

Downburst 

5.15 


" 

" 



13c 

" 

" 

" 

15.45 

15b 

" 

" 


15.45 

15c 

" 

" 


15.45 

17a 

Pilot 

S-Shape 

Longitudinal 

5.15 

17b 

" 

" 

" 

10.30 

17c 

" 

" 

■' 

15.45 

19b 

" 

" 

" 

15.45 

19c 

n 

" 

" 

15.45 

* 

Pilot 

1-Cosine 

Downburst 

5.15 



S-Shape 

Longitudinal 

5.15 

23a 

Pilot 

1-Cosine 

Downburst 

10.30 



S-Shape 

Longitudinal 

10.30 

23b 

' Pilot 

1-Cosine 

Downburst 

20.60 



S-Shape 

Longitudinal 

20.60 

23c 

Pilot 

1-Cosi ne 

Downburst 

20.60 



S-Shape 

Longitudinal 

20.60 


Case missing 

























X (in) 

a) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and frequency. 



b) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and 2 Wp^ frequency. 



X (m) 


c) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and 1/2 cop^ frequency. 

Figure 6 Computer-simulated approach trajectory of a B727 with fixed 
controls through a longitudinal sine wave wind. 
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(-38 sec) [3]. The encounter with the increasing head wind caused the 
aircraft to initially float above the glide slope. The encounter with 
the decreasing head wind to increasing tail wind portion of the wave 
caused a drastic decrease in airspeed of approximately 11 m/s (21.4 kts) 
causing the aircraft to drop below the glide slope. The aircraft then 
picked up airspeed during its descent below the glide slope and sub- 
sequently ascended above the glide slope. Thus, the encounter with the 
longitudinal shear wave excited an airspeed oscillation thereby causing 
flight path deviation. In this case, the airspeed ranged between 65 m/s 
and 88 m/s and angle of attack ranged between 3.7° and 6.1°. The air- 
speed deterioration parameters for this case were V+ = 5.02 m/s and V- = 
-5.93 m/s. Figure 6b shov/s the trajectory of the fixed-control B727 
encountering a 5.15 m/s (10 kt) head wind to tail wind sine wave at 2 cop^ 
(-19 sec). The aircraft remains in the wave for a shorter period of time 
than in the aforementioned case so that aircraft inertia is able to over- 
come some of the effect of the wave on airspeed. In this case, airspeed 
ranged between 67 m/s and 80 m/s and angle of attack ranged between 4.2° 
and 5.9°. The airspeed deterioration parameters were V+ = 2.78 m/s and 
V- = -3.68 m/s, which are of lower magnitude than for the phugoid period 
case. These parameters indicate lower airspeed variation for the 2 cjp^ 
wave than for the Wp^ wave, and therefore less flight path deviation. 
Figure 6c shows the trajectory of a B727 encountering a 5.15 m/s (10 kt) 
head wind to tail wind sine wave at 1 /2 03p^. The deviations from the 
glide slope are misleadingly large for this case because the aircraft 
remains out of trim for a longer period than for the two previous runs. 
However, the airspeed ranged between 71 m/s and 79 m/s and angle of 
attack varied from 4.4° to 5.2°. The airspeed deterioration parameters 
for this case were V+ = 1.9 m/s and V- = -2.37 m/s. Thus, the airspeed 
deterioration parameters were largest for the Wp^^ wave, indicating that 
the airspeed oscillations excited by the longitudinal shear wave will be 
largest for waves at the phugoid frequency. 

Figure 7 shows the computed trajectories with simulated pilot 
control through ajp^^, 2 cop^, and 1/2 tOp^^ sine waves, respectively. The 
ground speed/ vertical velocity control system I described in Section 2.3 
was used in these cases. It can be seen that the excessive flight path 
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a) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and o)p^ frequency. 



X (m) 


b) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and 2 frequency. 



X (m) 

c) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and 1/2 Wp^ frequency. 

Figure 7 Computer-simulated approach trajectory of a B727 using control 
system I through a longitudinal sine wave wind. 
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deviations of the preceding fixed-control runs have been dramatically 
lessened. In all cases, the aircraft remained on or above the glide 
slope for most of the run. A comparison between the cop^ and 2 copj^ 
cases shows less flight path and airspeed deviation for the 2 tOp^^ case 
than for the cop^^ case, which is to be expected following the results of 
the fixed-control trajectories. However, for the 1/2 Wp^ case, the air- 
craft ballooned above the glide slope and began to recover towards the 
end of the run. The ballooning and greater spread of airspeed deviation 
is due to the aircraft being in the wind for a longer duration, whereas 
the Wp^ and 2 cjp^ winds are relatively transient. 

4.3 Computed Aircraft Trajectory Results 

The following cases were run using the test plan described in 
Section 4.1. In all cases the B727 aircraft was trimmed (flaps 30°) for 
a 3.0° glide path angle with an approach airspeed of 70 m/s and an angle 
of attack of 6.2°. All wind profiles are encountered at x = 0.0 m. The 
flight profiles of the test pilot simulator runs have also been included 
in this section for comparison purposes. They are described in detail in 
Chapter 5. The four flight path deterioration parameters pertaining to 
each case are given on each trajectory plot. 

The airspeed/flight path angle control system II described in 
Section 2.3 was used in these cases. It should be noted that the fixed 
gains used in the formulation of this model represent an initial effort 
to model pilot control response to wind shear profiles and to simulate 
engine response characteristics. 

4.3.1 Sine Waves 

1. 5.15 m/s at Wp^^. Figure 8a shows the simulated trajectory of a 
B727 through a 5.15 m/s head wind to tail wind sine wave at Wp^. The 
aircraft balloons above the glide slope when encountering the increasing 
head wind portion of the wave. The pilot attempts to lower the thrust 
and pitch down to remain on the 3° glide slope. However, the aircraft 
next encounters the head wind to tail wind portion of the v/ave and is 
caught in a decreasing airspeed situation with low thrust and the nose 
pitched down. The pilot must then pitch the nose up and increase thrust 

26 



I 





a) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and frequency. 



b) Simulated trajectory through a longitudinal sine wave of 10.3 m/s 
amplitude and frequency. 



c) Simulated trajectory through a longitudinal sine wave of 15.45 m/s 
amplitude and frequency. 

Figure 8 Computer-simulated approach trajectory of a B727 using control 
system II through a longitudinal sine wave wind. 
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30B. 2 


GS+ = 1 .00 



X (m) 

a) Simulated trajectory through a longitudinal sine wave of 5.15 m/s 
amplitude and frequency. 



X (m) 


b) Simulated trajectory through a longitudinal sine wave of 10.3 m/s 
amplitude and Up^ frequency. 



c) Simulated trajectory through a longitudinal sine wave of 15.45 m/s 
amplitude and Wp^ frequency. 

Figure 9 Approach trajectory of a B727 aircraft performed by a test 

pilot in the NASA Ames simulator through a longitudinal sine 
wave wind. 
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to arrest the increased sink rate and to damp out the "roller coaster" 
effect of the longitudinal shear wave. The aircraft then reintercepts 
and drops below the glide slope at the termination of the run. However, 
it can be seen that the "roller coaster" effect is being damped out by 
pilot control. The amount of damping the pilot can apply is dependent 
upon the altitude at which the wave is encountered, which in this case 
is a relatively low 300 m. The airspeed varied from 63.7 m/s to 81.2 
m/s, angle of attack varied from 3.4° to 8.5°, and thrust ranged from 
idle to approximately 83 percent full forward. 

2. 10.3 m/s at 03^^^. Figure 8b shows the simulated trajectory of a 
B727 through a 10.3 m/s sine wave at Wp^. A comparison of the airspeed 
deterioration parameters reveals that the airspeed deviations are greater 
for the 10.3 m/s wave amplitude than for the 5.15 m/s amplitude, which 

is to be expected. The flight path deviation parameter spread [(GS+) - 
(GS-)] is larger for the 10.3 m/s wave than for the 5.15 m/s wave which 
is also to be expected. For this run, the airspeed ranged between 59.8 
m/s and 89.6 m/s and angle of attack varied between 1.8° and 10.0°. 

Thrust varied from idle to full forward. 

3. 15.45 m/s at Wp^. Figure 8c shows the flight path of a B727 
flown through a 15.45 m/s sine wave at Wp^^. The flight path deviation 
spread has again increased and the airspeed deviation parameters have 
also increased. Airspeed ranged between 56.9 m/s and 93.8 m/s and angle 
of attack varied between 1.1° and 11.4°. Thrust ranged from idle to full 
forward. 

4. 20.6 m/s at Wp^. The flight path of a B727 flown through a 

20.6 m/s longitudinal sine wave at Wpj^ is given in Figure 10a. Again, the 
GS+/GS- spread has increased and the airspeed parameters have increased. 
The airspeed varied between 54.4 m/s and 97.0 m/s, angle of attack ranged 
between 0.6° and 12.6° and thrust varied from idle to full forward. 

5. 20.6 m/s at 2 Wp^^. The trajectory of a B727 flown through a 
20.6 m/s sine wave at 2 Wp^ is given in Figure 10b. The aircraft under- 
went large airspeed deviations within a short time span and attempts to 
control the airspeed resulted in the aircraft entering the stall regime. 

Therefore, the trajectory is probably invalid past x = 1,000.0 m where 
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a) Simulated trajectory through a longitudinal sine wave of 20.6 m/s 
amplitude and frequency. 



X (m) 

b) Simulated trajectory through a longitudinal sine wave of 20.6 m/s 
amplitude and 2 frequency. 



X (m) 

c) Simulated trajectory through a longitudinal sine wave of 20.6 m/s 
amplitude and 1/2 cDp^ frequency. 

Figure 10 Computer-simulated approach trajectory of a B727 using control 
system II through a longitudinal sine wave wind. 
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X (m) 

a) Simulated trajectory through a longitudinal sine wave of 20.6 m/s 
amplitude and Wp^ frequency. 



b) Simulated trajectory through a longitudinal sine wave of 20.6 m/s 
amplitude and 2 Wp^ frequency. 



c) Simulated trajectory through a longitudinal sine wave of 20.6 m/s 
amplitude and 1/2 cop^ frequency. 

Figure 11 Approach trajectory of a B727 aircraft performed by a test 
pilot in the NASA Ames simulator through a longitudinal sine 
wave wind. 



the thrust was decreased to idle to bleed off the airspeed increase and 
the head wind to tail wind shear reached a maximum. Airspeed reaches a 
minimum of 45.4 m/s. This resulted in a descent rate and stall condition 
which probably would have been beyond the pilot's recovery ability due to 
the low altitude. Although the results of this case are contradictory to 
the hypothesis that the phugoid frequency waves are more severe than 
waves of other frequencies, the authors believe that a large wind gradi- 
ent due to large amplitude and short wave length may be detrimental to a 
safe approach. It appears that wave amplitude is an important factor in 
flight path deterioration within some frequency bandwidth bracketing the 
phugoid frequency. The extent of this bandwidth requires further inves- 
tigation. In summary, the large wind gradient and the attempt to over- 
come it clearly put the aircraft in a dangerous flight regime. 

6. 20.6 m/s at 1/2 Figure 10c shows the trajectory of a B727 

flown through a 20.6 m/s sine wave at 1/2 Wp^. The aircraft encounters 
the steadily increasing head wind and rises above the glide slope. The 
airspeed increase is soon arrested by idling the thrust even though the 
head wind is still increasing. Because of the amplitude, duration, and 
steadiness of the wave, the ground speed decreases from 70 m/s to 44 
m/s after 30 sec. As the aircraft enters the head wind to tail wind 
shear portion of the wave, the airspeed begins to drop and, as the air- 
craft drops below the glide slope, thrust is increased. However, the 
ground speed is also beginning to increase due to the tail wind and with 
the increased thrust, the pilot ends up overcontrolling towards the end 
of the run. Although this wind field clearly presents no problem with 
excessive airspeed oscillations and deviation below the glide slope, 
there is evidence of overcontrol by the control logic gains used to 
represent the pilot. The set response is too fast for such a long wave 
and the fixed gains for thrust and elevator control should be readjusted 
due to the quasi -steadiness and large wave amplitude of the wind field. 
The airspeed deviation parameters are small relative to the previously 
mentioned runs. The airspeed ranged between 63 m/s and 84 m/s, angle of 
attack varied between 2.9° and 8.7°, and thrust varied between idle and 
95 percent full forward. 
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In summary, the four flight path deterioration parameters showed 
increases with increased wave amplitude at the phugoid frequency. 

However, in the case of the 2 wave, the large wind gradient caused 
the aircraft to approach the stall regime. In the 1/2 Wp^ wave over- 
control was observed although it was noted that the airspeed oscillations 
were small for this case. 


4.3.2 1 - Cosine Downbursts 

1. 5.15 m/s at tOp^. The flight profile of a B727 flown through a 
5.15 m/s 1 - cosine downburst at Wp^ is shown in Figure 12a. As the air- 
craft encounters the down draft, it falls below the glide slope despite 
nose pitch-up and increased thrust. However, as with the longitudinal 
wave, a "roller coaster" oscillation was set up in vertical velocity 
resulting in deviations below the glide slope. It should be noted that 
the V+ and V- parameter values are lower than those of the 5.15 m/s 
longitudinal sine wave at Wp^. For this run, angle of attack remained 
fairly steady (5.0° to 7.4°), pitch varied between 3° and 7°, and air- 
speed ranged between 67 m/s and 76 m/s. 


2. 10.3 m/s at Wp^. The trajectory of the aircraft flown through 

a 10.3 m/s downburst appears in Figure 12b. Again, the "roller coaster" 
effect is noted and the aircraft is blown below the glide slope. The 
deterioration parameters are larger than those of the previous run. The 
angle of attack varied between 3.8° and 8.0°, pitch angle ranged from 
2.8° to 10.8°, and airspeed varied between 65.3 m/s and 78.7 m/s. Thrust 
setting reached idle and full forward limits. 


3. 15.45 m/s at ojp^^. In Figure 12c, it can be seen that the 15.45 

m/s downburst at the phugoid frequency forces the aircraft farthest below 
the glide slope. Pitch angle varied between 2.8° and 14.9° and angle of 
attack varied between 3.4° and 8.7°. Thrust setting ranged between idle 
and full forward settings and airspeed ranged between 63.0 m/s and 80.1 
m/s . 


4. 15.45 m/s at 2 co Figure 14b shows the flight profile of a 

B727 encountering a 15.45 m/s downburst wave at 2 Wp^. Because this wave 
is shorter in time duration than that of the phugoid period wave, the 
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a) Simulated trajectory through a vertical 1 - cosine wave of 5.15 m/s 
amplitude and frequency. 



b) Simulated trajectory through a vertical 1 - cosine wave of 10.3 m/s 
amplitude and oop^ frequency. 



c) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and Wp^ frequency. 

Figure 12 Computer-simulated approach trajectory of a B727 using control 
system II through a vertical 1 - cosine wave downburst. 
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X (m) 

a) Simulated trajectory through a vertical 1 - cosine wave of 5.15 m/s 
amplitude and frequency. 



X (m) 

b) Simulated trajectory through a vertical 1 - cosine wave of 10.3 m/s 
amplitude and Wp^ frequency. 



X (m) 

c) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and ujp^ frequency. 

Figure 13 Approach trajectory of a B727 aircraft performed by a test 
pilot in the NASA Ames simulator through a vertical 1 - 
cosine wave downburst. 
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a) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and frequency.* 



b) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and 2 frequency. 



c) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and 1/2 cOp^ frequency. 

Figure 14 Computer-simulated approach trajectory of a B727 using control 
system II through a vertical 1 - cosine wave downburst. 

*Repea ted from Figure 12 for consistency. 
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a) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and frequency.* 



X (m) 


b) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and 2 cjOp^^ frequency. 



X (m) 


c) Simulated trajectory through a vertical 1 - cosine wave of 15.45 m/s 
amplitude and 1/2 Wp^ frequency. 

Figure 15 Approach trajectory of a B727 aircraft performed by a test 
pilot in the NASA Ames simulator through a vertical 1 - 
cosine wave wind. 

*Repeat¥d from Figure 13 for consistency. 
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aircraft remains in the down draft for a shorter time and the altitude 
deviation from the glide slope is not as severe as it was in the phugoid 
wave. The airspeed deviation parameters are similar to those of the 
phugoid wave. The pitch angle varied from 2.7° to 16.0° and angle of 
attack ranged between 3.0° and 10.3°. Airspeed ranged from 60.5 m/s to 
81.8 m/s and idle to full throttle limits were reached. 

5. 15.45 m/s at 1/2 The flight of a B727 through a 15.45 m/s 

1 - cosine downburst at 1/2 tOp^ is shown in Figure 14c. The magnitude 
and duration of this wave forces the aircraft below the glide slope and 
causes premature impact. This is also due to the fact that the aircraft 
encounters the wave relatively close to the ground. Airspeed deviation 
is very small as indicated by the low values of V+ and V- and the fact 
that airspeed varied only between 66,6 m/s and 70.1 m/s. Angle of attack 
varied between 6.2° and 7.3° and pitch angle varied between 3.3° and 
14.6°. Thrust ranged from approach thrust to full forward. 

In summary, the flight profile deteriorates with increasing down- 
burst amplitude at the phugoid frequency. However, the worst case 
occurred not at the phugoid frequency but at the longest wave (1/2 Wp^). 
Thus it appears that the longer the aircraft is in the downburst, the 
greater the flight path deterioration. 

4.3.3 S-Shape Waves 

1. 5.15 m/s at Wp^. Figure 16a shows the simulated trajectory of a 
B727 through a 5.15 m/s longitudinal S-shape shear wave at Wp^^. This 
wave presents very little problem in flight path and airspeed deviations 
and is handled v/ell by the pilot. In this run, angle of attack varied 
between 5.6° and 6.9°, airspeed ranged between 68.3 m/s and 72.1 m/s, 
and thrust setting varied from 20 percent to 53 percent full power. 

2. 10.3 m/s at Wp^. The flight profile through a 10.3 m/s wave 
at the same frequency is presented in Figure 16b. Although the four 
deterioration parameters have increased, this wave does not present a 
serious control problem for the pilot. The angle of attack varied 
between 5.0° and 7.6°, airspeed deviated between 66.3 m/s and 74.4 m/s, 
and thrust ranged from idle to 73 percent full power. 
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a) Simulated trajectory through a longitudinal S-shape wave of 5.15 m/s 
amplitude and frequency. 



b) Simulated trajectory through a longitudinal S-shape wave of 10.3 m/s 
amplitude and Wp^ frequency. 



X (m) 

c) Simulated trajectory through a longitudinal S-shape wave of 15.45 m/s 
amplitude and Wp^ frequency. 

Figure 16 Computer-simulated approach trajectory of a B727 using control 
system II through a longitudinal S-shape wave wind. 
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a) Simulated trajectory through a longitudinal S-shape wave of 5.15 m/s 
amplitude and frequency. 



b) Simulated trajectory through a longitudinal S-shape wave of 10.3 m/s 
amplitude and cop^ frequency. 



X (m) 

c) Simulated trajectory through a longitudinal S-shape wave of 15.45 m/s 
amplitude and Wp^ frequency. 

Figure 17 Approach trajectory of a B727 aircraft performed by a test 

pilot in the NASA Ames simulator through a longitudinal S-shape 
wave wind. 
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3. 15.45 m/s at Wp^. In the case of a 15.45 m/s wave at Wp^ 

(Figure 16c), the deterioration parameters have increased and control is 
more difficult to attain. Angle of attack was maintained between 4.4° 
and 8.0° and airspeed ranged between 65.3 m/s and 76.4 m/s. Thrust 
setting varied from idle to 87 percent full forward. 


4. 15.45 m/s at 2 o)p^. Figure 18b shows the flight path of the air- 

craft flown through a 15.45 m/s S-shape wave at 2 Wp^. It appears that 
the larger longitudinal wind gradient (due to large wave amplitude and 
short wavelength) and the pilot's control inputs cause greater airspeed 
oscillations and a greater drop below the glide slope. In this case, 
angle of attack varied from 3.1° to 9,0° and airspeed ranged between 
62.1 m/s and 81.8 m/s. Idle and full thrust limits were reached. 


5. 15.45 m/s at 1/2 Wp^. The aircraft encounter with the 15.45 m/s 

wave at 1/2 ojp^ caused little airspeed and flight path deviation. The 
slight overcontrol tendencies of the pilot control system logic due to 
the fixed gains and the long wave duration can be noted near the end of 
the trajectory (Figure 18c). The airspeed oscillated between 69.2 m/s 
and 71.4 m/s and angle of attack was maintained between 5.8° and 6.5°. 
Thrust ranged between 27 percent and 55 percent full forward. 

In summary, the flight profiles deteriorated with increasing wave 
amplitude at the phugoid frequency. However, the effects of the' large 
wind gradient caused the largest deterioration in the case of the 2 ojp^^ 
wave at 15.45 m/s amplitude. The 1/2 ojp^ was seen to cause no serious 
control problems. 

4.3.4 Com bi nation S-Shape Waves and Downbursts 

The previously mentioned longitudinal S-shape and sine waves and 
vertical downburst profiles were examined individually to determine the 
aircraft/pilot response to that particular wind component. The next four 
cases were run to examine the combined effects of longitudinal S-shape 
waves and 1 - cosine downbursts. This combined profile is the most 
realistic with respect to the actual thunderstorm downburst cell wind 
field. Fujita and Caracena's analysis [2] of the Eastern 66 accident in 
Figure 20 shows the downburst cell located along the aircraft's approach 
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a) Simulated trajectory through a longitudinal S-shape wave of 15.45 m/s 
amplitude and frequency.* 



b) Simulated trajectory through a longitudinal S-shape wave of 15,45 m/s 
amplitude and Z o)p^ frequency. 



c) Simulated trajectory through a longitudinal S-shape wave of 15.45 m/s 
amplitude and 1/2 Wp^ frequency. 

Figure 18 Computer-simulated approach trajectory of a B727 using control 
system II through a longitudinal S-shape wave wind. 

*Repeated from Figure 16 for consistency. 
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X (m) 

a) Simulated trajectory through a longitudinal S-Shape wave of 15.45 m/s 
amplitude and frequency.* 



b) Simulated trajectory through a longitudinal S-shape wave of 15.45 m/s 
amplitude and 2 frequency. 



c) Simulated trajectory through a longitudinal S-shape wave of 15.45 m/s 
amplitude and 1/2 Wp^ frequency. 

Figure 19 Approach trajectory of a B727 aircraft performed by a test 
pilot in the NASA Ames simulator through a longitudinal 
S-shape wave wind. 

*Repeated from Figure 17 for consistency. 
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-pi 


M FT 

213.4 ( 700 ) 

182.9 ( 600 ) 

152.4 ( 500 ) 

121.9 ( 400 ) 

91.4 ( 300 ) 
61 .0 ( 200 ) 

30.5 ( 100 ) 

0.0 ( 0 ) 

Figure 20 The flight path of Eastern 66 and associated wind field as 
analyzed by Fujita and Caracena [2]. 




I 

path. As the downburst reaches the ground, it spreads out. The aircraft 
encounters a head wind to tail wind shear causing airspeed loss and a 
strong down draft causing descent below the glide slope. Thus, the com- 
bined effects of the longitudinal and vertical winds in the thunderstorm 
downburst cell present a serious hazard to an aircraft operating at low 
altitude. 

1. 5.15 m/s at A comparison between Figure 12a (5.15 m/s 

downburst at and Figure 21 (combination of downburst and S-shape 

wave) reveal little difference in flight path deterioration. Angle of 
attack varied between 4.6° and 7.6° and airspeed varied between 66.2 m/s 
and 75.4 m/s. Thrust setting ranged between idle and 75 percent full 
forward power. 

2. 10.3 m/s at Comparison of Figure 12b (10.3 m/s downburst 

at ojp^) and Figure 22a (combination of downburst and S-shape wave) also 
shows little difference in flight path deterioration. This reveals that 
the smaller amplitude 5.15 m/s and 10.3 m/s S-shape waves, which were 
shown to have little effect by themselves, cause no additional difficulty 
when acting in conjunction with the smaller amplitude downburst waves. 

3. 20.6 m/s at Wp^^. The combined airspeed loss due to a 20.6 m/s 
S-shape longitudinal wave and vertical velocity increase due to a 20.6 m/s 

1 - cosine downburst at w . forced the aircraft to descend below the 

ph 

glide slope and impact prematurely (Figure 22b). The large value of V- 
indicates a large drop in airspeed below the approach airspeed. The 
airspeed in this case decreased to a minimum of 54.3 m/s and the angle 
of attack reached 13.2°. 

4. 20.6 m/s at 2 Wp^^. Figure 22c shows the trajectory for the 

2 o)pj^ case. The large wind gradients and combined effects of the S- 
shape wave and downburst caused the aircraft to descend below the glide 
slope. In this case, the airspeed dropped to 49.8 m/s and the angle of 
attack entered the stall regime (approximately 16°) resulting in a loss 
of control. The large value of airspeed warning parameter V- indicates 
a large airspeed drop below the approach value and would warn the pilot 
of the possibility of a stall. 
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X (m) 

Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 - cosine wave of 5.15 m/s amplitude and frequency. 


Figure 21 Computer-simulated approach trajectory of a B727 using control 
system II through a combination longitudinal S-shape and 
vertical 1 - cosine wave wind. 



X (m) 

a) Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 - cosine wave of 10.3 m/s amplitude and frequency. 



X (m) 

b) Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 - cosine wave of 20.6 m/s amplitude and Wp^ frequency. 



X (m) 

c) Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 - cosine wave of 20.6 m/s amplitude and 2 Wp^ frequency. 

Figure 22 Computer-simulated approach trajectory of a B727 using control 

system II through a combination longitudinal S-shape and vertical 
1 - cosine wave wind. 
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a) Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 ~ cosine wave of 10.3 m/s amplitude and frequency. 



b) Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 - cosine wave of 20,6 m/s amplitude and Wp^ frequency. 



c) Simulated trajectory through a combination longitudinal S-shape and 
vertical 1 - cosine wave of 20.6 m/s amplitude and 2 frequency. 

Figure 23 Approach trajectory of a B727 aircraft performed by a test 
pilot in the NASA Ames simulator through a combination 
longitudinal S-shape and vertical 1 - cosine wave wind. 
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Finally, in Figure 24 a test run was made through a 15.45 m/s sine 
wave at using the pilot control system logic III as discussed in 
Section 2.3. This system controls airspeed by throttle rate control and 
pitch attitude and touchdown point by elevator control . Comparison with 
the flight path deterioration parameters in Figure 8c reveals that this 
pilot control model does not allow the aircraft to sink as far below the 
glide slope. Additionally, the airspeed deterioration parameters V+ and 
V- are reduced. However, a relatively high vertical velocity is needed 
to bring the aircraft to the touchdown point due to a rise above the 
glide slope at a relatively low altitude. 
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X (m) 


Figure 24 Computer-simulated approach trajectory of a B727 using control 
system III through a longitudinal sine wave wind of 15.45 m/s 
amplitude and to ^ frequency. 


CHAPTER 5. DESCRIPTION OF FLIGHT SIMULATOR TESTS 
CONDUCTED AT NASA AMES 


5.1 Flight Simulator Results 

The following tests corresponding to the test plan described in 
Table 1 were run in the B727 flight simulator at NASA Ames. The profiles 
were flown by a test pilot for the purpose of evaluating the response of 
the aircraft/pilot system to various thunderstorm environment wind shears. 


5.1.1 Sine Waves 


and 11 


1. 5.15 m/s at This run did not cause serious control diffi- 

culty for the pilot. The pilot decreased thrust and pitched nose down 
when encountering the increasing head wind portion of the wave. As the 
head wind sheared to a tail wind, the pilot pitched nose up and increased 
thrust. The pilot maintained better airspeed control (67.0 m/s to 74.2 
m/s) as compared with 63.7 m/s to 81.2 m/s in the computer simulation. 
Flight path angle and angle of attack were also kept within closer con- 
straints. Thrust control ranged from idle to 88 percent power compared 
to idle and 83 percent in the computer model. 


2. 10.3 m/s at Wp^. This simulation presented a greater control 
problem than the previous run. The flight path and pitch angle varia- 
tions were larger than those of the previous run and angle of attack 
ranged from 2,5° to 10.5° as compared with 1.8° and 10.0° for the compu- 
ter simulation. Airspeed varied from 58.8 m/s to 77.3 m/s (59.8 m/s to 
89.6 m/s for computer model) and thrust ranged from idle to full as in 
the computer run. 

3. 15.45 m/s at Wp^. This simulation was similar in control diffi- 
culty to the 10.3 m/s run. This is probably a result of the pilot learn- 
ing from the preceding simulations. The airspeed was maintained between 
61.3 m/s and 85.0 m/s as compared to 56.9 m/s and 93.8 m/s in the compu- 
ter run. Near idle to full forward thrust limits were reached. 


51 



4. 20.6 m/s at od^. 


In this run, the pilot showed further evidence 


of "learning" since the wave frequency had been kept the same as for the 
previous runs. However, greater thrust and elevator control were needed 
to maintain the flight path. The pilot kept airspeed between 61.9 m/s 
and 89.6 m/s compared with 54.4 m/s and 97.0 m/s in the computer simula- 
tion. Angle of attack reached 8.2° compared with 12.6° in the computer 
simulation. 


5. 20.6 m/s at 2 This profile took the pilot by surprise and, 

as in the computer simulation, aircraft control was lost. However, the 
test pilot was able to recover before landing but not without significant 
altitude and flight path angle deviations. Angle of attack reached 18.5° 
(stall regime) as in the computer model. Airspeed fell to 50.0 m/s 
(45.4 m/s in simulation). This simulation agrees well with the computer 
model and suggests that large wind gradients due to large wave amplitude 
and short wavelength may seriously affect aircraft operations. The 
effects of large amplitude at frequencies within some bandwidth of the 
phugoid frequency must be further examined. 

6. 20.6 m/s at 1/2 This wave presented no difficulty with 

regard to flight path angle and deviation below the glide slope. Air- 
speed was kept between 64.9 m/s and 85.6 m/s (63.0 m/s and 84.0 m/s for 
computer run). However, tighter pitch and angle of attack control by 
the test pilot than by the simulated pilot made the difference in more 
accurate flight path control. It was noted that the computer- simulated 
pilot tended to overcontrol near the termination of the run and was 
therefore not an accurate representation of pilot performance in a large 
amplitude, long wavelength shear. 

5.1.2 1 - Cosine Downbursts (Figures 13 and 15) 

1. 5.15 m/s at cop^. Flight through this profile caused slight 

flight path angle deviations. The pilot pitched nose up to a maximum of 
7° (7° in computer simulation) and increased thrust during the downburst 
encounter. Angle of attack remained fairly constant between 4.2° and 
6.4° compared with the model's 5.0° and 7.4°. Airspeed varied from 67.0 
m/s to 71.1 m/s compared to 67.0 m/s to 74.0 m/s for the model. 
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2. 10.3 m/s at cop^. Larger pitch and thrust were needed to control 
the flight path angle in this profile. Pitch ranged between 2.5° to 
11.0° (compared with the simulation's 2.8° to 10.8°) and maximum angle of 
attack reached 6.3° compared with 8.0° for the model. Airspeed varied 
from 66.5 m/s to 70.0 m/s (model, 65.3 m/s to 78.7 m/s). The test pilot 
performed better on this run than the simulated pilot. 

3. 15.45 m/s at Wp^. Large flight path angle and significant alti- 
tude deviations occurred during this run. The pilot had to make large 
pitch and thrust corrections to readjust and maintain the approach. A 
maximum pitch angle of 13.0° was reached (14.9° for computer model) and 
maximum angle of attack reached 7.4° (8.7° for computer simulation). 
However, airspeed and thrust were more accurately controlled by the test 
pilot than by the simulated pilot and he was therefore able to recover 
where the computer model could not. 

4. 15.45 m/s at 2 Wp^. As in the computer simulation, the test 
pilot has less difficulty with this downburst wave than with the 1/2 Wp^^ 
wave. Large flight path and altitude deviations occurred over the 19 sec 
duration of the wave and large thrust and pitch inputs were needed to 
restabilize the approach. Maximum angle of attack reached 9.3° (10.3° 
for model) and airspeed varied from 64.4 m/s to 74.7 m/s (model, 60.5 m/s 
to 81,8 m/s) . 

5. 15.45 m/s at 1/2 cup^. As in the computer simulation, the test 
pilot had the greatest difficulty with this profile. The test pilot was 
better able to negotiate this downburst but not without large thrust and 
pitch activity (14.5° maximum, 14.6° for model). The test pilot was able 
to keep airspeed between 67.0 m/s and 70.1 m/s and angle of attack nearly 
constant at 5.0°. However, large flight path angle and altitude devia- 
tions occurred before the aircraft was restabilized. 

5.1.3 S-Shape Waves (Figures 17 and 19) 

1. 5.15 m/s at Wp^. The pilot encountered little difficulty with 

this profile. The computer model also experienced little difficulty. 
Airspeed varied from 67.0 m/s to 70.1 m/s compared with 68.3 m/s and 
72.1 m/s for the model. 
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2. 10.3 m/s at Wp^. This profile also presented no problem for the 

test pilots As airspeed decreased due to the head wind to tail wind 
shear, thrust was increased and then decreased when the aircraft was 
retrimmed for the remaining tail wind. Airspeed was held between 67.0 
m/s and 70.6 m/s compared with 66.3 m/s and 74.4 m/s for the model. The 
test pilot flew better than the model pilot because of better airspeed/ 
thrust management. 


3. 15.45 m/s at (a .. The pilot had little problem with this 

S-shape wave while the simulated pilot experienced larger airspeed and 
altitude deviations. Despite the fact that the test pilot has thus far 
flown more precisely than the model pilot, the precision of this approach 
may also be due to "learning." Airspeed ranged between 66.5 m/s and 70.6 
m/s compared with 65.3 m/s and 76.4 m/s for the model. Flight path, 
pitch, and angle of attack were kept nearly constant. 


4. 15.45 m/s at 2 This shear wave required a large thrust 

input over a short time period by the pilot. The maximum angle of attack 
was 8.2° (9.0° for the model) and airspeed ranged between 62.4 m/s and 
72.0 m/s compared with 62.1 m/s and 81.8 m/s for the model. Although the 
test pilot flew this approach better than the model pilot due to better 
thrust/airspeed control, this case presented the greatest control problem 
for both. This supports the reasoning that large wave amplitudes at 
frequencies within some bandwidth bracketing the phugoid frequency may 
cause hazardous airspeed and altitude deviations. 


5. 15.45 m/s at 1/2 Wp^. This profile required very little thrust 

and pitch input from the test pilot as with the -simulated pilot. Air- 
speed varied from 67.5 m/s to 70.4 m/s compared with 69.2 m/s and 71.4 
m/s for the model. Flight path, pitch, and angle of attack were closely 
maintained by the test pilot and the model, although the test pilot again 
flew a more precise approach. 

5.1.4 Combination S-Shape Waves and Dov/nbursts (Figure 23) 


1. 5.15 m/s at Wp^. The test pilot encountered little difficulty 

with this profile. When the downburst and longitudinal shear were 
encountered, the pilot increased pitch attitude and thrust to maintain 
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approach. This profile did not cause very significant problems for the 
computer simulation although the test pilot flew it more precisely. 

2. 10.3 m/s at This profile required larger pitch (12.5®) and 

thrust changes by the test pilot. Some altitude and flight path angle 
deviations occurred and a relatively safe approach was made again with 
greater precision than that of the simulated pilot. 

3. 20.6 m/s at Wp^. As in the computer model, this profile caused 
very large altitude and flight path angle deviations. The test pilot 
stated that the simulated aircraft had been "forced significantly below 
the glide slope" before he was able to barely recover. Airspeed reached 
a minimum of 57.4 m/s compared with 54.3 m/s for the model and maximum 
angle of attack reached 13.2° (13.2° for the model). 

4. 20.6 m/s at 2 Wp^. This profile caused loss of control as did 
the previous case. Again significant flight path angle, airspeed, and 
altitude deviations occurred resulting in a loss of control. This was 
shown to occur in the computer simulation as well. Airspeed reached a 
minimum of 50.8 m/s (49.8 m/s for the model) and the angle of attack 
reached a maximum of 19.1°. However, the test pilot eventually managed 
to bring the aircraft under control. 

5.2 Comparisons Between Computer Results 
and Simulator Results 

The trajectory program was tested for sine wave and S-shape longitu- 
dinal shear waves, 1 - cosine downbursts and combinations of longitudinal 
and vertical waves. The results of these simulations were compared to 
simulations through the same wind fields flown in the NASA Ames simulator 
by a test pilot. In terms of aircraft/pilot response, the computer model 
compared well with the simulator for the sine waves, S-shape waves, down- 
bursts, and combinations. However, some discrepancies existed with 
regard to the degree of flight path and airspeed control between the 
computer model and the test pilot. Although the control logic for the 
model pilot is similar to the control strategy of the test pilot, the 
test pilot flew consistently better than the model pilot. This is due 
to the fixed gain structure of the computer model pilot. 
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The fixed gains for elevator and thrust control were selected as 
reasonable control responses to airspeed, flight path angle, and pitch 
attitude variations. A real pilot does not behave in such a rigid 
manner. In reality, a pilot acts in a variable gain decision-making 
process which is probably not adequately included in the simplified 
model. This fixed gain structure of the model allows for lower pilot 
damping of the flight path and airspeed oscillations induced by encounter 
with wave disturbances. The test pilot is clearly of better skill than 
the computer model pilot. In addition, the test pilot had the opportu- 
nity to "learn" the types of profiles he was flying during the tests. 

The airspeed deterioration parameters calculated for the fixed- 
stick cases presented earlier in Figure 6 showed that the V+ and V- 
values for the small amplitude wave (5,15 m/s) at the phugoid frequency 
were larger than those of the 2 and 1/2 Wp^ waves at the same ampli- 
tude. A small amplitude wave was used for the illustration since the 
aircraft could not handle waves of larger amplitudes in the fixed-stick 
mode. However, in the cases of large amplitude waves and pilot control, 
the computer and simulator runs showed that the large amplitude 15.45 
m/s S-shape and 20.6 m/s sine waves at 2 Wp^ presented the greatest 
control difficulty for the aircraft/pilot system. The airspeed deterio- 
ration parameters V+ and V- calculated from the computer simulations 
(Figures 25 and 26) are shown to increase with increasing longitudinal 
wave amplitude for waves at the phugoid frequency. The largest V- 
values were attained for the 2 Wp^ waves. However, the V+ value of the 
20.6 m/s sine wave at 2 Wp^ is not as large as that of the phugoid 

frequency wave, which is misleading. The values of V+ and V- from the 

computer simulations were smallest for the 1/2 tOp^ S-shape and sine 
waves. In comparison with the computer results, similar trends can be 
noted from the airspeed deterioration parameters from the simulator runs 
for the S-shape waves and for V+ values of the sine waves. However, in 

the case of the sine waves, the simulator values of V- tend to be 

misleading due to the pilot "learning" the profiles and "fine tuning" 
his control procedures. 

Control difficulty was encountered by the computer model and test 
pilot in flight through downbursts (particularly the 1/2 cop^ wave) and 
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a) V+ values versus sine wave amplitude and frequency. 



b) V- values versus sine wave amplitude and frequency. 


Figure 25 Airspeed deviation parameters for longitudinal sine wave winds. 
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a) V+ values versus S-shape wave amplitude and frequency. 



a) V- values versus S-shape wave amplitude and frequency. 


Figure 26 


Airspeed deviation parameters for longitudinal S-shape wave 
winds . 


58 


in the combination S-shape longitudinal waves and 1 - cosine downbursts. 
For the downbursts, the computer-generated and simulator airspeed dete- 
rioration parameters (Figure 27) were lowest for the 1/2 wave which 
caused the largest flight path deviation. This is reasonable since the 
long wave downburst does not have a pronounced effect on airspeed 
deviation but instead causes the aircraft to descend below the glide 
slope with the steadily descending air mass. The aircraft remains in 
this long wave for 76 seconds. Therefore, airspeed deterioration is 
probably not a meaningful warning parameter for application to downbursts. 
It may be noted that the glide slope deviation parameter GS- (Tables 2 
and 3) for downbursts show low values corresponding to large descent 
below the glide slope for the 1/2 waves and that this type of warning 
parameter might be more meaningful. 

For the combination S-shape longitudinal waves and 1 - cosine down- 
bursts, the decreasing airspeed and the descending air mass forcing the 
aircraft below the glide slope presented the most difficulty for the 
simulated pilot and test pilot. The V+ values given in Figure 28a show a 
misleading trend. However, Figure 28b shows that the largest V- values 
for airspeed deviation correspond to the worst control cases and would 
therefore provide a meaningful warning of hazardous wind conditions. 
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a) V+ values versus 1 - cosine wave amplitude and frequency. 



b) V- values versus 1 - cosine wave amplitude and frequency. 


Figure 27 Airspeed deviation parameters for vertical 1 - cosine wave 
winds. 
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a) V+ values versus combination S-shape and 1 - cosine wave amplitude 
and frequency. 



(m/s) 

(m/s) 


b) V- values versus combination S-shape and 1 - cosine wave amplitude 
and frequency. 


Figure 28 Airspeed deviation parameters for combination longitudinal 
S-shape and vertical 1 - cosine wave winds. 
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CHAPTER 6. CONCLUDING REMARKS 


The feasibility of a system for warning and detection of longitudi- 
nal wind shear has been discussed. The aircraft trajectory program, 
pilot control models, and the flight path deterioration parameters have 
also been discussed. The results of the computer and simulator runs 
through the longitudinal S-shape waves and sine waves indicate a need 
for further studies to examine the effects of large longitudinal wind 
gradients due to large amplitude waves at short wavelengths. A param- 
etric analysis of a broader range of wave frequencies must be carried 
out to determine the bandwidth which is most hazardous to aircraft 
operations. 

Control difficulties were noted with 1 - cosine downbursts, particu- 
larly the long duration wave at large amplitude and the strong downbursts 
combined with longitudinal shear. The authors believe that the combina- 
tion longitudinal S-shape and vertical 1 - cosine wind profile is the 
most realistic model of a downburst cell wind field in the vicinity of 
the ground. Therefore, additional simulator studies should be carried 
out to study the response of the aircraft/pilot system to a wider range 
of combination downbursts and longitudinal shear waves. 

The authors believe that the Doppler radar is the optimum wind mea- 
surement device for a wind shear detection system which provides the 
aircraft with advanced warning of hazardous shear conditions. As men- 
tioned above, a downburst cell wind profile consists of a vertical 1 - 
cosine wave component and a longitudinal S-shape wave component. Because 
of the interrelationship of the longitudinal and vertical winds, it is 
probable that the measurement of only the longitudinal wind component 
along the flight path and the detection of longitudinal shear waves may 
provide a sufficient warning when a downburst hazard exists. 

The present computer model has the capability to simulate pilots of 
varying skill. Further studies can be performed on "fine tuning" the 
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pilot model to simulate the best, average, or worst pilot. A simulator 
study of pilot elevator and thrust inputs in flight through wind gradi- 
ents of varying magnitudes would yield necessary data for use in the 
simplified numerical model. It is also possible to put in the actual 
autopilot used in the aircraft to serve as the baseline pilot. The 
autopilot characteristics could easily be entered into the control logic 
in the aircraft trajectory model and used as the standard pilot for the 
wind shear warning and detection system. 

In addition, the aircraft trajectory model used in the actual 
warning system would require a more accurate but equally simplified 
model of engine response to pilot throttle input. The most accurate and 
concise aerodynamic data for commercial transports must also be available 
for improved accuracy of the warning system aircraft trajectory software. 

Finally, the issue of how flight path deterioration should be quan- 
titatively determined and presented as a go/no-go signal to the pilot or 
air traffic controller must be resolved. The results of the computer 
study shows that while the altitude deviation parameters can be mislead- 
ing (except possibly for downburst warnings), the airspeed deviation 
parameters may yield useful results when longitudinal disturbances are 
present. The authors believe that V- may be a useful warning parameter 
because it reveals information about how far below an aircraft may 
deviate. An example of the usefulness of V- as a warning parameter is 
given by the two severe combination S-shape and downburst waves which 
took the test pilot by surprise, causing serious control difficulty. 

The large values of V- would have warned the pilot that the aircraft 
airspeed would deviate significantly below and into the stall 
region. This would indicate a no-go situation and the warning system 
would have alerted the aircraft before entering the Doppler-monitored 
flight path. Further tests are needed to determine the go/no-go decision 
value. It is also possible to use the minimum airspeed obtained during 
the simulated trajectory as a quantitative deterioration parameter. If 
this value was below the go/no-go decision value (some value around the 
stall speed) the aircraft would be forewarned of the hazardous conditions 
existing along its anticipated flight path. 
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APPENDICES 



APPENDIX A 


LANDING CHARACTERISTICS OF A B727 CLASS MEDIUM-SIZED JET TRANSPORT* 

Flaps 30° 

Gear Down 

Glide Slope Angle -3.0° 

V trim airspeed, 70.0 m/s 

m aircraft mass, 63,958.0 kg 

6 2 

lyy moment of inertia, 6.1 x 10 kg - m 

c mean chord length, 4.57 m 

2 

A wing area, 145.0 m 

Cj_ 0.74 

^0 

Cj_ 6.99/rad 

a 

C, 0.361/rad 

^E 

Cj_ 10.0/rad 

q 

C^^ -7.6/rad 

a 

Cp 0.152 

^0 

Cp 0.3/rad 

Cn”* 2.4/rad^ 

U 2 

C|> 0.0/rad 

C -0.25 

">0 

C -1.40/rad 

m 

a 

C -1.59/ rad 

m p 

C -30.0/rad 

C -2.16/rad 

m* 


*Aerodynamic data used in an airline flight simulator. 
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APPENDIX B 


SAMPLE CALCULATION OF FLIGHT PATH DETERIORATION PARAMETERS 


The flight path deterioration parameters (FPDP) quantitatively 
describe the aircraft trajectory by glide slope deviations and airspeed 
deviations from reference airspeed. Normalized altitude, HP/HG, was 
chosen as the FPDP representing glide slope deviation, where HP is the 
height of the aircraft above the ground and HG is the height of the 
glide slope above the ground. This formulation was chosen over root 
mean square calculations of glide slope deviation (aH = HP - HG) 
because it accounts for the fact that altitude deviations become more 
serious at lower altitudes than at higher altitudes. 

The airspeed deviation parameter chosen was simply the deviation of 
airspeed from the reference value (aV = V - Vf)- I" both glide slope 
and airspeed deviation parameters, the positive deviations are examined 
separately from the negative deviations to avoid the cancellation errors. 

Figure B.l shows a hypothetical aircraft trajectory of which a quan- 
titative analysis of glide slope and airspeed deviation is desired. The 
values of airspeed and altitude for each point on the trajectory is 
given in Table B.l. Also given in Table B.l are the FPDP's, airspeed 
deviation from reference airspeed, and normalized altitude. Using the 
data from Table B.l, the FPDP's are calculated as follows: 


Computation of V+ 


Total points airspeed is above or at = 15 = j 


V + 


1 j 
= T ^ 

J i=l 


(V^ - 


Vref) = 


1 

15 


0.0 + 1.0 + 2.0 + 4.0 + 5.0 + 6.0 + 4.0 
+ 1.0 + 0.0 + 1.0 + 2.0 + 3.0 + 1.0 
+ 0.0 + 0.0 


^ [30.0 m/s] 


2.0 m/s 


deviation above 


(average airspeed 
= +2.0 m/s) 
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Figure B.l Hypothetical aircraft trajectory. 


(T> 

<£> 


Point 

Airspeed 
V (m/s) 

1 

70.0 


71.0 


72.0 


74.0 


75.0 


76.0 


74.0 

8 

71.0 

9 

69.0 

10 

68.0 

11 

67.0 

12 

69.0 

13 

70.0 

14 

71.0 

15 

72.0 

16 

73.0 

17 

71.0 

18 

70.0 

19 

70.0 

20 

69,0 

Total 20 



TABLE B.l Hypothetical Aircraft Trajectory Data 
























Computation of V- 


Total points airspeed is below = 5 = k 

" F i, ■ ''ref^ " F ( 

^ [-8.0 m/s] = I -1 .6 m/s~i (average airspeed 
deviation below = -1.6 m/s) 

Computation of GS+ 

Total points aircraft is above or on glide slope = 14 = n 


1 1 
■ n HGT " 14 


' 1.00 + 1.05 + 1.06 + 1.08 + 1.07 + 1.05 + 1.05 1 
+ 1.03 + 1.00 + 1.07 + 1.20 + 1.43 + 1.20 
+ 1.50 


[15.79] = 1 1 . 13| (average displacement of aircraft 
above glide slope = 13 percent above glide slope 
height) 


Computation of G$- 


Total points aircraft is below glide slope = 6 = m 


m HP. 


GS- = I 
m - ^ 


1 


i=l 


HG. 


1 

6 


0.97 + 0.93 + 0.86 + 0.84 + 0.86 + 0.89 


0.89 


(average displacement of aircraft below glide 
slope = 11 percent below glide slope height) 
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